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Abstract 



A simple analysis technique is described which allows to intercalibrate the response 
of imaging atmospheric Cherenkov telescopes in stereoscopic telescope arrays at a 
^ ; level of 1-2%. 

m 

^ ■ With the next-generation imaging atmospheric Cherenkov telescopes currently 

O . under construction, stereoscopic arrays of telescopes such as VERITAS [1], 

^ \ CANGAROO III [2] and H.E.S.S. [3] play a dominant role. Stereoscopic sys- 

Qh| tems of Cherenkov telescopes for VHE gamma-ray astrophysics provide - com- 

P ' pared to a single telescope - superior angular resolution, energy resolution and 
background rejection, as demonstrated by the HEGRA group [4,5]. Making 

^ . use of information concerning the impact point of showers and the height of 



the shower maximum, energy resolutions around 10% are possible [6]. To ac- 
tually obtain this resolution, the relative response of the telescopes, i.e. the 



' relation between the light yield and the digital counts provided by the ADC 

^ system must be known with very good precision. Photon detectors within one 

telescope can be fiat-fielded using a light pulser in the center of the dish. 
Telescope-to-telescope intercalibration is more complicated, and recently e.g. 
the use of cosmic ray trigger rates and the so-called throughput method has 
been advocated for this purpose [7]. In the following, a simple technique for 
telescope intercalibration is discussed (see also [8]). 

Once a telescope is properly fiat-fielded, one needs only one global calibra- 
tion factor per telescope, relating the image size - the sum of the amplitudes 
of all image pixels - to the photon yield, up to some overall calibration fac- 
tor common to all telescopes of an array. The response of two telescopes can 
be compared by selecting showers with cores in the middle between two tele- 
scopes, and comparing the image sizes of both telescopes. However, in systems 
with more than two telescopes, one needs to be careful not to introduce a bias 
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due to other telescopes involved in the triggering or the event reconstruc- 
tion. Testing the technique with the HEGRA telescope system, the following 
procedure was followed: 



• To compare the response of two telescopes, events were selected where both 
telescopes under consideration had triggered (and possible others in addi- 
tion) 

• Only information from the two telescopes under consideration was used in 
reconstruction the shower; images from other telescopes were ignored. 

• To reduce the influence of trigger thresholds, a cut well above the threshold 
was applied; the sum of the image sizes of the two telescopes had to be 
above 200 photoelectrons. 

• To improve image quality and reconstruction, only images within a radius 
of 1.75° from the camera center were used, avoiding image truncation. The 
(stereo) angle between the two image axes had to be at least 30°. Recon- 
structed shower directions had to be within 0.7° from the telescope axis. 



The relative response of two telescopes i and j was then characterized by 
plotting the size asymmetry Ug versus the asymmetry in impact distance 
(Fig. 1), where 
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and 

ar{ij) 



Ti + Ti 



Here, Si is the image size in telescope i and n its distance to the recon- 
structed shower impact point, measured perpendicular to the shower axis. 
The asymmetry parameters were chosen because they treat both telescopes in 
a symmetrical fashion; one could also use log(sj/sj) or similar quantities. To 
evaluate the relative response of the two telescopes, one needs to determine 
the average Ug for Ur — 0. This can be achieved by averaging Ug for bins of 
and then fitting a smooth curve to < as(ar) >, or, if event statistics is high 
enough, by simply cutting on a^. For the following examples, the mean value 
of as{ar ~ 0) were obtained by fitting a Gaussian to the distribution of for 
\ar\ < 0.05. The method was first applied to Monte-Garlo simulations of the 
HEGRA telescope system in its 1997 configuration, with four telescopes; one 
of the outer telescopes was still incomplete. In the simulation, the response of 
the four telescopes (for historical reasons labelled CT3, CT4, CT5, CT6) was 
adjusted as 1 : 0.8 : 1.1 : 1.3. The "measured" size asymmetries are listed in 
the following table; the last column gives the values expected on the basis of 
the input response factors. 
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Fig. 1. Distribution of size asymmetry versus distance asymmetry a^, for a tele- 
scope pair of the Mkn 501 data set. The dashed hnes indicate the cut range used 
to derive 0^(0^ ~ 0). 



Telescopes 


asymmetry 
±0.002 


true value 


CT3-CT4 


0.113 


0.111 


CT3-CT5 


-0.052 


-0.048 


CT3-CT6 


-0.128 


-0.130 


CT4-CT5 


-0.160 


-0.158 


CT4-CT6 


-0.242 


-0.238 


CT5-CT6 


-0.082 


-0.082 



Prom the first three hnes, one readily determines a relative response of the 
telescopes to 1 : 0.797 : 1.109 : 1.292 ^ , in agreement with the input values 
to better than 1%. With six measurements for three independent quantities 
(arbitrarily defining the response factor for CT3 as unity), one should deter- 
mine optimum calibration factors by a fit to all measurements; such a fit yields 
response ratios of 1 : 0.797 : 1.104 : 1.300, with typical errors of 0.003 and 
a of 5.2 for 3 degrees of freedom, indicating that small systematic effects 



^ Table and results rounded to three digits; calculations and fits used higher 
precision 
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may be present beyond the statistical errors. 

Tiie metiiod was then apphed to actual HEGRA data, from observations 
of Mkn 501 in 1997. Two independent data samples were considered, (a) a 
background-subtracted gamma-ray sample obtained by cutting on the mean 
scaled width of images and on the reconstructed direction relative to the 
source, and (b) a cosmic-ray sample. The measured asymmetries are given 
in the following table: 



Telescopes 


sample (a) 


sample (b) 




±0.006 


±0.003 


CT3-CT4 


0.033 


0.021 


CT3-CT5 


0.071 


0.073 


CT3-CT6 


0.117 


0.114 


CT4-CT5 


0.037 


0.039 


CT4-CT6 


0.089 


0.084 


CT5-CT6 


0.023 


0.045 



Fits for the telescope response factors yield 

sample (a): CT3 : CT4 : CT5 : CT6 = 1 : 0.937 : 0.857 : 0.797 

sample (b): CT3 : CT4 : CT5 : CT6 = 1 : 0.948 : 0.872 : 0.798 

with typical errors of 0.01 for sample (a) and 0.005 for sample (b). The results 
are consistent within errors; the values for the fits are 9.7 and 7.5, again 
slightly worse than expected for purely statistical errors. 

In summary, the method allows an easy intercalibration of imaging atmo- 
spheric Cherenkov telescopes in a telescope system. Both the simulations and 
the comparison of two independent data samples shows that a precision of 
1-2% can be reached, which is ample even for a 10% overall energy resolution. 
For systems with more than two telescopes, the equation system is overdeter- 
mined and consistency checks are possible, e.g. in the form of the value of 
an overall fit. The method can be applied to systems with arbitrary numbers 
of telescopes, provided that they do not separate into multiple distant clusters 
of telescopes. The absolute calibration for the whole system can be determined 
in a subsequent step, e.g., using cosmic-ray detection rates (see, e.g. [7,9]). 
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